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Electron-phonon interaction at the Be(0001) surface. 
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We present a first principle study of the electron-phonon (e-p) interaction at the Be(0001) surface. 
The real and imaginary part of the e-p self energy (E) are calculated for the F surface state in the 
binding energy range from the F point to the Fermi level. Our calculation shows an overall good 
agreement with several photoemission data measured at high and low temperatures. Additionally, 
we show that the energy derivative of ReT, presents a strong temperature and energy variation close 
to Ef, making it difficult to measure its value just at Ef- 
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For many years now a considerable effort has been 
directed toward the study of electron-phonon (e-p) ef- 
fects and quasiparticle dynamics in low dimensional 
systems. Important examples include semiconductor 
heterostructures, superconducting materials with lay- 
ered structures and electrons in two-dimensional sur- 
face states. Experimentally, e-p effects are becom- 
ing increasingly accessible on surfaces, through powerful 
techniques like high resolution angular-resolved photoe- 
mission spectroscopy (PES), scanning tunnelling spec- 
troscopy (STS), and time-resolved two-photon photoe- 
mission spectroscopy (TR-2PPE). 

Be(0001) offers in this respect a textbook study case. 
Here the surface states form a fully metallic, high density 




km r 




f km r 



FIG. 1: Top panel shows the electronic band structure of 
a fully relaxed 12 layer Be(0001) slab, the F surface state is 
indicated by the green line. Bottom panel shows the phonon 
dispersions where the red line represents the Rayleigh mode. 
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FIG. 2: Different contributions to the Eliashberg function for 
the hole state a) at Ef, b). at F point. 



two-dimensional electron gas, rather well decoupled from 
the bulk Be substrate, itself only of semi-metal character 
0, 0] • The surface states increase the density of electron 
states (DOS) at the Fermi energy (Ep) by a factor of 4-5 
in the surface layer . This radical change of character 
between surface and bulk goes along with some other very 
special properties of the Be(OOOI) surface 0- Following 
the DOS arguments, the surface A parameter, representa- 
tive of the e-p interaction, is believed to be about a factor 
of 4-5 larger than in bulk Be 0,H,IS0- Values of A rang- 
ing from 0.65 to 1.2 were obtained from the high temper- 
ature surface state broadening 0, or by fitting the real 
part of the low temperatures surface state self-energy to 
simple models ya, la, |]j . The scatter between these values 
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FIG. 3: a) The e-p contribution to i?eS (orange) and 7mE° 
(black) at T — as function of binding energy. Jm£° is 
broken up into contributions from the intraband scattering 
(green) and contribution from the Rayleigh mode (red), b) 
Black line represents the total value of Ao (Eq. |3J and green 
and red lines the contributions from intraband scattering and 
Rayleigh mode, respectively. 



is disturbing. Also, a value as large as 1.2 may raise ques- 
tions concerning the stability of the normal metallic state 
of this surface. On the other hand, it is well known that 
simple model calculations, such as Debye models, are 
surely too crude, neglecting band structure effects and 
the important contribution of the surface phonon modes 
0. Electron energy-loss spectroscopy (EELS) measure- 
ments showed that the surface phonon dispersions at 
Be(OOOl) do qualitatively differ from those obtained by 
using truncated bulk models . 

In this letter we present the first parameter free ab- 
initio calculation of the e-p interaction at a real metal 
surface, Be(OOOl), that contains all the ingredients en- 
tering the e-p interaction, namely, surface electron and 
phonon states, as well as electronic screening. Our de- 
tailed analysis results in excellent agreement with all ex- 
perimental data. We also motivate more future experi- 
mental work. 

A fully relaxed 12 layer slab has been considered. We 
have used the density-functional theory in the local den- 
sity approximation (LDA). The e-p matrix elements and 
phonon frequencies were calculated using density func- 
tional perturbation theory 0], through the PWSCF pro- 
gram 11] . A norm conserving pseudopotential was used 
with a cutoff of 22Ry for the expansion of the electronic 
wave function in plane waves. For the Brillouin zone elec- 
tronic structure integrations an order one Hermite-Gauss 
smearing technique was used, with smearing parameter 
a = 0.05 Ry and as many as 102 special k points in the 
irreducible wedge of the 2D Brillouin zone (IWBZ). The 



linear response calculation for phonon modes 0] have 
been performed with 30 special q points in the IWBZ, 
yielding surface phonons in excellent agreement with ex- 
periments [jj and previous calculations [l^ . 

The basic quantity in the study of the e-p interaction is 
the so called Eliashberg function (ELF) This func- 
tion measures the phonon density of states weighted by 
the e-p interaction. It represents the probability of emit- 
ting a phonon of energy uj at T = 
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The ELF contains all the information needed for the cal- 
culation of the real and imaginary part of the self-energy. 
In Eq. □ g i jff + =(^ i n \AscFV^\^ f M? ) is the e-p matrix 
element where 'J are the electronic state wave functions 
and AscfV^ denotes the change of the self-consistent 
electron potential induced by a phonon mode with mo- 
mentum q and mode index v\ i refers to the electronic 
state (surface state) index, / to the final state, E denotes 
the binding energy, Ni(E) the i state partial density of 
states at energy E and Hbz is the area of the 2D Bril- 
louin zone. 

The imaginary part of the self-energy related to the 
lifetime of the excited electronic state through t _1 = 
2 ImYi can be written in terms of the ELF, 



/>oo 

ImE°(E;T) = n du a 2 F l {uj,E) 
Jo 

(l-f{E-u)+f(E + u) + 2n(u)), 



(2) 



where / and n are the Fermi and Bose distribution func- 
tions, respectively. The real part reflects renormaliza- 
tion of the electronic band due to the interaction with 
phonons and can also be calculated from the ELF 13], 

/OC /"OC 
dv I duj a 2 F t (uj,E) 
-00 Jo 
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For the calculation of the ELF, the k integration was 
carried out with 102 special points in the IWBZ and with 
30 points for the q integration. The two Dirac deltas 
containing the electron band energies were again replaced 
by order one Hermite-Gauss smearing functions. In order 
to gain confidence in the convergence of the calculation, 
we performed a set of calculations with different smearing 
parameters in the range 0.0025 — O.lRy, which gives an 
estimated error of ~ 5 % for the surface state broadening 
at all energies except inside a window of 0.4 eV around 
Ep where we estimate an error of ~ 15 %. 

Fig. n snows the electron and phonon band structures. 
In the upper panel of Fig. ^ the f point surface state 
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FIG. 4: a) Surface state imaginary part of the self-energy 
near Ef- ImT, corresponds to the unrenormalised result 
(Fig. the renormalised JmE is represented for temper- 
atures 12K (violet) and 40K (orange). The inset shows 
modified dispersions at 12K and 40K. b) ReT, near Ef 
at 12K and 40K (orange). The inset shows the quantity 
A(E,T) = — 9 ■ Note that at low temperature (12K) 

A(E F ,T) ~ 0.9 is consistent with Ao(£f) ^ 0.9 from Fig. El 



band is sticked out (green line), near the M point another 
surface electronic state also appears. In the lower panel, 
the phonon dispersions are shown. The Rayleigh wave 
(red line) is a surface localised phonon mode polarised 
mainly along the surface normal. This and other surface 
phonon modes are visible within the surface band gaps. 

Fig. [3shows the calculated ELF for a hole at Ep (up- 
per panel) and at the f point (lower panel). The ELF 
is broken up into contributions corresponding to the in- 
traband scattering (green line) and the Rayleigh mode 
scattering contribution (red line). From the upper panel 
of Fig. [3 the decay of the hole at Ep is dominated by 
the intraband scattering contribution (green line). The 
sum of the contributions from the decay into M point 
surface state (not represented) and the intraband contri- 
bution (green line) gives ~ 90% of the total (solid line) 
at Ep. Thus, the decay process of the hole at Ep is 
sharply confined in space to almost exclusively the sur- 
face state region. By contrast, the character of the inter- 
action changes radically for the hole at the f point (lower 
panel), for which the decay into bulk states dominates. 

Photoemission experiments enable the measurement 
of the high temperature behaviour of the broadening, 
T(E) — 2nkBTXn(E), for a state with arbitrary bind- 
ing energy E 0, . Theoretically A (E) is determined 
by the first reciprocal moment of the ELF 01 • 
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Eq. ^indicates an enhanced role of the ELF at low ener- 
gies. Fig. [2] shows that the weight of the surface phonon 
modes is concentrated around lod/2 (lud — 80 meV), 
thus a fitting procedure of Ao parameter through a De- 
bye model (a 2 F = X (lu/ujd) 2 ) will give rise in this case 
to an excessive low temperature broadening compared to 
a more realistic calculation using the same Debye energy 
Ma and Ao. The top panel of Fig. shows the real (or- 
ange) and imaginary part (black) of the self-energy cal- 
culated for T=0 and as function of the binding energy. 
The imaginary part is resolved into contributions coming 
from intraband scattering (green) and the contribution of 
the Rayleigh mode (red) . The imaginary part varies by a 
factor of two going from Ep to the T point. The bottom 
panel of Fig. shows the variation of the Ao parameter 
defined through Eq. 0] as function of the binding energy. 
From this figure we deduce that the temperature deriva- 
tive of the high-T surface state broadening is twice for 
the hole at the Ep compared to that at the T point. 

Recently, the identity between the partial derivative of 
i?e£ at Fermi level and X (Ep), that is A = X (Ep) — 
A(E F ,T = 0), 



(5) 



has been used to determine the e-p coupling parameter 
A [MIEH • However, it is of crucial importance here to 
distinguish A and Ao- The results of Fig. [3] do not take 
into account the renormalization of electron bands due to 
the real part of the self-energy. We denote by E% the pure 
LDA unrenormalised bands and by E^ the renormalised 
ones including i?e£, obtained by solving numerically the 
equation E k = E° k + ReT,(E k ) (see Fig. 0J. It is well 
known that close to Ep the electron velocity, e-p coupling 
form factor and the lifetime itself are renormalised • 
Provided that the system is isotropic and \ A(E, T)\ « 
1 , an estimation of the renormalised imaginary part [l3| 
can be obtained from 
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where ImS denotes the renormalised quantity. While the 
breakdown of the approximation expressed by Eq. is 
clear for values of A(E, T) close to -1, this form generally 
helps us to understand how the deformation of the elec- 
tronic band caused by the real self-energy affects itself 
the lifetime of the surface state. 

In the top panel of Fig. 21 we show the unrenormalised 
imaginary part of the self-energy (7m£°, dotted line) cal- 
culated from Eq. [21 and the renormalised ones at T = 40 
(orange) and at T = 12K (violet). The result for 12-K" is 
displayed only for reasonable values of 1/(1 + A). These 
two temperatures, 12K and 40if, are the same as those 
considered experimentally, Refs. Q and 0] , respectively. 
In the bottom panel, the real part of the self-energy at 
T = 12K (violet) and at T = 40K (orange) are depicted. 
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The inset of the bottom panel shows the derivatives of 
the real part of the self-energy at these two temperatures. 

In Ref. |{| a detailed measurement of ImT and ReT 
is reported for T = 4QK. An 'unusual' peak of ImT 
for the hole at binding energies around the maximum 
phonon energy (~ u>d) was observed. This behaviour 
is very well described by our calculation at the same 
temperature (orange). The calculated real part of the 
self-energy also shows excellent agreement with the mea- 
surement, and our A(Ep, T = 40AT) ~ 0.7, is in excellent 
agreement with their measured 0.7 ± 0.1 (Inset of the 
bottom panel of Fig. 0J. In Ref. 0, Q the derivative 
of the real self-energy part was presented for T — 12 K 
suggesting a value of A(E F ,T = 12K) ~ 1.15. As ex- 
plained in Ref. this value was slightly overestimated 
due to the fixed angle photocmission experiment, and a 
corrected value of 0.87 was suggested instead. Our cal- 
culated A(Ep, T — 12K) ~ 0.9 is in excellent agreement 
with this corrected value (inset of the bottom panel of 
Fig. 0J. Ref. |j] reported high temperature measure- 
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FIG. 5: Temperature dependence of the broadening at Wd/2, 
at ujd and at 3/2ujd binding energies. Black lines represent 
the renormalised broadening according to Eq. ©. Note the 
minimum at temperatures ksT ~ 13 meV. Black dotted line 
for ujd indicates region where Eq E]is not reliable. The inset 
shows the ImT, dependence on binding energy (Eb). 

ments of the broadening for the surface state hole at bid- 
ing energy ~ 350 meV, obtaining an estimate Ao(E=350 
meV) of 1.15 ±0.1 , which compares reasonably well with 
a calculated Ao(E= 350meV) ~ 0.9 (see bottom panel of 
Fig. El. 

Analysing the temperature dependence of the broad- 
ening in more detail, we note that for larger T , the max- 
imum value of the T-derivative of ReT decreases, and 
so does the effect of the renormalising denominator in 
Eq. 6 (Fig. 0}. By contrast, increasing temperature 
favours phonon occupation: thus we have two compet- 
ing tendencies and a minimum can be expected in the 
temperature dependence of ImT,. Detailed calculation of 
ImT predicts a minimum of ImT around T ~ 150A' for 
states close to ujp, binding energy (see Fig. [SJ. It would 
be interesting to verify this minimum in the tempera- 
ture dependence of the broadening. This effect should 
be even more pronounced for binding energies close to 



~ ujd where the effect of the renormalisation is stronger. 
Here in fact the group velocity selectively rises above its 
bare value, and so does the low temperature hole decay 
rate, giving rise to a local maximum as a function of hole 
energy. Removal of this effect by temperature causes the 
local minimum predicted in Fig 5. 

In summary we presented the first parameter-free 
ab — initio calculation of the e-p interaction in a surface 
state on a real metal surface, Be(0001). The surface state 
renormalisation and its broadening were calculated, and 
the latter was found to agree well with existing exper- 
imental data. The energy derivative of ReT presents a 
strong energy and temperature dependence near Ep , so 
that an accurate measurement of A = A{Ep) = Xo(Ep) 
requires temperatures T < 12A and energies very close 
to Ep. Measurements done under different conditions 
can produce very different A values which explains the 
variety of values in literature. Additionally, the effect of 
the renormalisation of the surface state band on its own 
broadening is pointed out and new experiments are sug- 
gested to confirm a new minimum predicted in the tem- 
perature dependence. Other conceivable consequences 
of a medium size surface state e-p coupling, including 
the possible stabilisation of a 2D BCS state at T=0 |l4j 
and/or the factors hindering it, remain to be explored. 

We acknowledge financial support from the Basque 
Government, the Max Planck Research Award Funds, 
the Spanish Ministerio de Education y Ciencia (MEC), 
and the Basque Country University. Work in SISSA was 
supported through MIUR COFIN and FIRB projects, by 
INFM/F and by the 'Iniziativa Trasversale Calcolo Par- 
allelo' of INFM. 



[1] R. A. Bartynski et al, Phys. Rev. B 32, 1921 (1985). 

[2] E. V. Chulkov et al, Surface Sci. 188, 287 (1987). 

[3] E. W. Plummer and J. B. Hannon, Prog. Surf. Sci. 46, 
149 (1994); L. I. Johansson, et al, Phys. Rev. Lett. 71, 
2453 (1993); P. T. Sprunger et al, Science 275, 1764 

(1997) . 

[4] T. Balasubramanian et al, Phys. Rev. B 57, R6866 

(1998) . 

[5] S. LaShell et al, Phys. Rev. B 61, 2371 (2000). 
[6] M. Hengsberger et al, Phys. Rev. Lett. 83 592 (1999). 
[7] M. Hengsberger et al, Phys. Rev. B 60, 10796 (1999). 
[8] A. Eiguren et al, Phys. Rev. Lett. 88, 066805 (2002). 
[9] J. B. Hannon et al, Phys. Rev. B 53, 2090 (1996). 
[10] S. Baroni et al, Rev. Mod. Phys. 73, 515 (2001). 
[11] S. Baroni et al, http://www.pwscf.org 
[12] M. Lazzeri et al, Surf. Sci. 402 715 (1998). 
[13] G. Grimvall, in The Electron- Phonon Interaction in Met- 
als, Selected Topics in Solid State Physics, edited by E. 
Wohlfarth (North-Holland, New York, 1981). 
[14] E. Tosatti, in Electronic surface and interface states on 
metallic systems, eds, E. Bertel and M. Donath (World 
Scientific, Singapore, 1995), p. 67. 



